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ABSTRACT 



X 

We survey all possible supersymmetric three-body decays of the top quark in 
the framework of the MSSM and present detailed numerical analyses of the 
most relevant cases. Although the two-body channels are generally dominant, 
it is not inconceivable that some or all of our most favourite two-body SUSY 
candidates could be suppressed. In this event there is still the possibility that 
some of the available three-body SUSY modes might exhibit a substantial 
branching fraction and/or carry exotic signatures that would facilitate their 
identification. Furthermore, in view of the projected inclusive measurement 
of the top-quark width T t in future colliders, one should have at one's dis- 
posal the full second order correction (electroweak and strong) to the value of 
that parameter in the MSSM. Our analysis confirms that some supersymme- 
tric three-body decays could be relevant and thus contribute to Ft at a level 
comparable to the largest one loop supersymmetric effects on two-body modes 
recently computed in the literature. 



1. Introduction 



The era of the top quark has just begun |]. To a large extent we were prepared to 
enter the long-announced new epoch and in the meanwhile a tremendous amount of work 
has piled up on top-quark observables[} Yet the extremely rich potential phenomenology 
and the far-reaching consequences that top-quark dynamics may have on the final picture 
that will emerge of the Standard Model (SM) of the electroweak and strong interactions 
definitely demands a new fully fledged wave of theoretical and experimental endeavor in 
Particle Physics. 

The SM has been a most succesful framework to describe the phenomenology of the 
strong and electroweak interactions for the last thirty years Q . The top quark itself stood, 
at a purely theoretical level -namely, on the grounds of requiring internal consistency, 
such as gauge invariance and renormalizability- as a firm prediction of the SM since the 
very confirmation of the existence of the bottom quark and the measurement of its weak 
isospin quantum numbers |J . Nevertheless, despite all the successes there are still too 
many questions left unanswered by the SM, especially in connection with the nature of 
the spontaneous symmetry breaking mechanism (SSB) and the purported existence of the 
extremely elusive Higgs particle, whether realized as a truly elementary particle or as an 
effective (composite) field. Due to the huge mass of the top quark, one expects it to be 
the most preferential fermion to which the Higgs particle should couple. Therefore, top 
physics is deemed to be the ideal environment for Higgs phenomenology. 

Lately the SM has been exposed to an escalate of experimental information potentially 
challenging some of its predictions to an unprecedented level. We are referring to the long- 
standing conundrum involving the high precision Z-boson observable Rb and some related 
observables . The issue about Rt> is specially offending, for it seems to consolidate with 
time-the present day discrepancy with the SM being at the 3.4 a level 0. Whether this 
anomaly is linked to an incomplete understanding of the experimental analysis of Z decays 
into b-quarks, or it can be licitly associated with some sort of physical effects beyond the 
SM, has not yet been established and at present there is a lot of controversy about it @ . Be 
as it may, from the theoretical point of view one is tempted to believe that physics might 
be taking a definite trend beyond the SM. One possibility is to look at the predictions of 
the supersymmetric (SUSY) extension of the SM. In this paper we take our chances in 
favour of the elementary Higgs particle interpretation of the SSB and we adhere to the 
supersymmetric extension of the SM, more specifically to the Minimal Supersymmetric 
Standard Model (MSSM) ||. In particular, there is in the literature quite a big amount 
of early || as well as of very recent SUSY work on Ri, ||10||-||14j|. In some of these works 

1 See e.g. Refs. [0, || and references therein. 
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it is shown that the discrepancies, although they cannot be fully accounted for, at least 
they can be significantly weakened under suitable conditions ]ll], [12], [L3| . 

In view of the new wave of SUSY potentialities, it is natural to study all possible phe- 
nomenological consequences that may be expected on supersymmetric top-quark physics; 
after all, the interactions between the top-quark and Higgs sector are doubled in a SUSY 
framework and one may hope that top-quark physics can be a window to both Higgs 
bosons and supersymmetric particles. As a first step in this direction one would like to 
assess the importance of real and virtual SUSY effects on top-quark decays. Here, too, a 
fairly respectable amount of work is scattered over the literature]]: 

i) Supersymmetric two-body decays of the top quark have been described at the tree 
level e.g. in Refs.g], 0-0; 

ii) Supersymmetric Higgs corrections to the conventional top-quark decay mode t — > 



W + b have been computed in Ref. p2 



iii) Supersymmetric electroweak quantum effects on t — > W + b mediated by the roster 
of genuine supersymmetric particles, such as sleptons, squarks, charginos and neutralinos, 
have been accounted for in Refs. [f23|, p4]j ; 

iv) Supersymmetric QCD corrections to t — > W + b have also been studied in detail in 
Ref. [H; 

v) Supersymmetric QCD effects on the charged Higgs decay of the top quark, t — > 
H + b, are generally relevant and can be spectacularly large in favourable regions of the 
MSSM parameter space [^, ^7|1 ; 

vi) The full plethora of electroweak quantum effects on the unconventional mode t — > 
H + b in the MSSM are also recently available and can be rather significant |T6| , |28| . 

We see that on the theoretical side there is a large amount of work ready to be used 
by experimentalists. Now, what about the prospects for an experimental measurement of 
these effects?. On the one hand, the measurement of Y(t — > W + b) will be carried out at 
the Tevatron at a level of ~ 10% and will be further reduced at the LHC In this respect 
we remind that the top-quark phenomenology is expected to significantly improve at the 
Tevatron 0, ^9] on the basis of a projected ten- fold increase of the luminosity via the 
Main Injector and Recycler facilities, together with a ~ 35% increase of the production 
cross-section at a 2TeV running energy (Run II), as compared to the fruitful 1.8 TeV 
past run (Run I). However, in a hadron machine one aims more at a measurement of 
specific top-quark production vertices, which are obviously related to the corresponding 
top-quark partial decay widths. For instance, one of the main goals at the Tevatron for 



Run II is the measurement of the single top-quark production cross-section [£9], [30fl , which 
gives essential information on the SM vertex tbW and, a fortiori, on the value of the 



2 For a review of some of these results, see Refs. fill 16 
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CKM matrix element Vtb- In the absence of new physics this measurement is equivalent 
to a determination of the top quark width. However, in the presence of new interactions 
beyond the SM, such as e.g. SUSY interactions, one may expect significant changes in 
the prediction for the cross-section which can be related to the hadronic partial widths of 
the Higgs bosons of the MSSM [[H], ^, 

On the other hand, from the point of view of an inclusive mo del- independent measure- 
ment of the total top-quark width, T t , the future e + e~ supercollider should be a better 
suited machine. In an inclusive measurement, all possible non-SM effects would appear on 
top of the corresponding SM effects already computed in the literature |34j] . As shown in 



Ref . |3"5| , one expects to be able to measure the top-quark width in e + e~ supercolliders at 
a level of ~ 4% on the basis of a detailed analysis of both the top momentum distribution 
and the resonance contributions to the forward-backward asymmetry in the tt threshold 
region. 

Clearly, for a consistent treatment of the corrections to the observable T t at second 
order of perturbation theory (in the strong and electroweak gauge couplings) one should 
include the tree-level contributions from all possible three-body decays of the top quark 
in the MSSM. As it happens, the contribution of some of these three-body decays turn 
out to be comparable to the largest SUSY quantum effects on the two-body channels 
mentioned above. Furthermore, it could occur that our most cherished SUSY two-body 
decays are not kinematically allowed or are significantly suppressed in some regions of 
parameter space. Therefore, supersymmetric three-body decays could be relevant and a 
detailed study is in order. Such a study is, to our knowledge, missing in the literature and 
it is precisely the task that we have undertaken in this article. The paper is organized 
as follows. In Section 2 we present an overview of two-body and three-body decays of 
the top quark in the MSSM. The Lagrangian interactions relevant to these decays in the 
MSSM are given in Section 3. The numerical analysis of the various partial widths, with 
special emphasis on the dominant modes, is detailed in Section 4. Finally, Section 5 is 
devoted to a discussion of the results, as well as of the possible signatures for the favourite 
decays, and we deliver our conclusions. An appendix is provided at the end to display 
some lengthy formulae. 

2. Decays of the top quark in the MSSM 

The weighted average of the CDF and DO determinations of the top-quark mass reads 



as follows 



m t = 175 ±9 GeV. (1) 

Due to the large mass of the top quark, there is plenty of phase space available for two- 
body and multibody decays. Within the minimal SM, the leading standard decay of the 
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top quark is the "canonical" decay 



t^W + b. (2) 
The associated partial width is given at the tree-level by 

ra.r^ri) = f-°0 1M a'/*(i, 4, ^% 

\8nV2j m t mf mf 

x \M 2 w {m\ + m 2 b ) + (m 2 - m 2 b ) 2 - 2M* W ] , (3) 

where 

A 1/2 (l, x 2 ,y 2 ) = ^[l-( x + y y][l-( x - y y]. (4) 

Additional standard decays to other quarks are of course possible but, contrary to the 
canonical decay (in which Vtb — 0.999, for three generations), they are CKM-suppressed. 
Numerically, for m t = 175 GeV (and m;, = 5 GeV) one has T(t — > W + b) ~ 1.54 GeV, 
which is very big as compared to Aq C d ~ 300 MeV, and therefore the top quark decay 
is basically a free-quark decay |37| . At this point we should recall that the latest mea- 
surements of the canonical branching ratio at the Tevatron still give room enough for top 



quark decays beyond the SM, namely they could reach BR{t — > "new") ~ 40% |36 



As for the SUSY two-body decays PI , jl7fl- Plfl i the leading modes are the following. 
On the SUSY-QCD side the top quark can disintegrate, if there is phase space enough, 
into the lightest stop (ti) and gluino (g), 

t^hg, (5) 

and on the SUSY electroweak side we have 

t^hxt, (6) 

and 

t^iiXi, (7) 

where b\ stands for the lightest sbottom and xf (xT) f° r the lightest chargino (neutralino). 
In some cases the next-to-lightest neutralinos can also enter and be relevant (see Sections 
4-5). Another conspicuous electroweak decay of the top quark in the MSSM involves a 
super symmetric charged Higgs in the final state, 

t^H + b. (8) 

This decay, if kinematically allowed, could be very promising and it has recently been 



studied up to one- loop order in great detail in Ref . p8| . Numerically, all the two-body 



SUSY modes can be important as compared to the canonical mode (0) in certain regions 
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of the MSSM parameter space. We shall compare their contribution with that of the 
leading three-body decays in Section 4. 

Concerning the three-body decays, there are the SM modes where the W + is a virtual 
particle that subsequently goes into lepton or quark final states. However, since the W + 
can be real, the two-body mode (0) followed by the decay of the W + as a real particle 
is overwhelming in the SM. In fact, this is the way in which the top quark has been 
discovered |l|] . Therefore, we are only interested in three-body decays of the top quark 
in the MSSM other than the three-body SM decays. We shall christen these decays, the 
SUSY three-body decays of the top quark. In contrast to the SM case, in the MSSM 
not all of them need to be suppressed as compared to the two-body modes, as we shall 
see. In an extreme situation, the three-body decays could be the leading SUSY decays 
of the top quark. There are a fairly big number of them, but in the end only a few can 
be of some relevance. In the following, and unless stated otherwise, we shall impose the 
following condition: the relevant three-body decays under study are only those decays in 
which the virtual particle is heavy enough that the corresponding SUSY two-body decay 
is kinematically forbidden. Later on we shall relax this condition in some especial cases. 

The following decays are essentially ruled out by our conditions or by direct phe- 
nomenological constraints: 

1. t — > iiW + X\'- Forbidden by phase space: 

m~ h + M x - < m t -M w ~ 100 GeV 

mr+M v - > 120 GeV. ^ 

where the second relation reflects the recent LEP 130/140 bound on the lightest 



chargino and sfermion masses [|38 



m h ,M x ± > 60 GeV, (10) 



2. t — > bix\W + : In this case we require 



m~ bi + M x o <m t - M w ~ 100 GeV 

m~ tl + M x o > m t (11) 



so that 



m- bi + M x + > m t 



M x o < 40 GeV 

M v + > 80 GeV + M Y o , 

Xi Xi ' 



(12) 



which cannot be fulfilled in the gaugino-higgsino window (/i, M) -see Section 3- of 
the MSSM parameter space, as we have checked. 
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3. t — > t\ H ± Xi : Impossible in the MSSM where m H + > M w , and so we cannot fulfil 
the phase space constraint 

m H+ <m t - m k - M x - < 60 GeV , (13) 

not even in the presence of SUSY radiative corrections, which would lower (only 
slightly) the charged Higgs boson mass P§f . 

4. t — > H + X\bi, t — > H + gh\: For the first process, we have the condition 

m H + <m t — M x o - m~ bl , (14) 

which forces the charged Higgs mass to be near the lowest limit allowed in the 
MSSM. For the second, we may admit of both light gluinos m- g = Oil) GeV |4(| or 
heavy gluinos [|4[] , 

m- g > 100 GeV. (15) 

In the light gluino case, since 

rriH+ < m t — m g — m~ hx < 120 GeV , (16) 

the two-body decay t — > H + b would already be allowed. In the heavy gluino sce- 
nario, the decay at stake is excluded since it would enforce an unacceptably light 
charged Higgs: 

m H + < m t — m g — ^ 20 GeV . (17) 

5. t — > tiW + Xi- It would require 

m h + M x - <m t -M w ~ 100 GeV , (18) 
which is ruled out by eq. (p~Q|) . 



Among the decays in the complete list of SUSY three-body decays of the top-quark 
in the MSSM which cannot be discarded by trivial arguments, we remark the following: 

I. t^br + v T 

II. t-^h°bW + 

III. t^b a W + g 

IV. t^bxlxt 

V. t^t b b a b 

VI. t—>b a cs b 

VII. t^b a r + V T 

VIII. t^bgxt 

IX. t^bt b b a 

X. t-ftf+iv, (19) 



where for each decay we sum over all sparticle indices allowed by phase space (see Section 
4). The decays quoted above do not exhaust the list of potential three-body modes, 
but the related modes not included in the list are less significant. For instance, decay I 
above proceeds both via a virtual W + and via a virtual charged Higgs boson, H + . The 
reason why this decay has been singled out over similar decays involving light quarks, 
e.g. t — > bud, is because for the latters the charged Higgs couplings to light quarks are 
suppressed as compared to the coupling to a r-lepton. As another example, consider decay 
II. It involves the lightest CP-even SUSY Higgs boson, h° ||17|| . Clearly, two related modes 
are obtained by replacing h° with the heavy CP-even Higgs boson, H°, or with the CP- 
odd Higgs boson, A . In the first case, the additional decay will obviously be suppressed 
by phase space; and in the second case, since the charged Higgs H + is prescribed to be 
heavy enough in order to prevent the two-body mode t — > H + b from occurring, it follows 
from the usual Higgs mass relations in the MSSM IL7]] that the A must be heavier than 



h°, and therefore the corresponding decay is phase-space obstructed. 

The Feynman diagrams contributing to the various decays ( |T9| ) are given in Figs. 1-8. 
Each process has been thoroughly studied and the main numerical results are provided 
in Section 4. The upshot of our analysis is that there are a few selected decays in the 
list ( |19D that could be of interest. As for the remaining decays, even though they cannot 
be dismissed by trivial arguments of the sort used in the cases considered before, they 
eventually prove to be irrelevant. 

3. Lagrangian interactions for top quark decays in 
the MSSM 

Although the Lagrangian of the MSSM is well-known H, it is always useful to project 
explicitly the relevant pieces and to cast them in a most suitable form for specific purposes. 
Even with this arrangement, some care is to be exercised in actual calculations, due to the 
Majorana nature of the neutralinos and the complicated mixing among the various fields. 
We shall perform our calculations in a mass-eigenstate basis. One goes from the weak- 
eigenstate basis to the mass-eigenstate basis via appropriate unitary transformations. 
Two classes of SUSY particles enter our computations: the fermionic partners of the 
weak-eigenstate gauge bosons and Higgs bosons (called gauginos, B, W, and higgsinos, 
H, respectively) and on the other hand the scalar partners of quarks and leptons (called 
squarks, q, and sleptons, I, respectively, or sfermions, /, generically) . Among the gauginos 
we also have the strongly interacting gluinos, g, which are the fermionic partners of the 
gluons. Within the context of the MSSM, we need two Higgs superfield doublets with 
weak hypercharges Yip = Tl- The (neutral components of the) corresponding scalar 
Higgs doublets give mass to the down (up) -like quarks through the VEV < >= v\ 



8 



^< #0 >= v 2 ). The ratio 



tan 5=—, (20) 



is a most relevant parameter throughout our analysis. 

Due to the "crural" L-R mixing between weak-eigenstate sfermions of a given flavor, 

q' a = Wi = ql, q'2 = qr} > (2i) 

there is a matrix relation with the corresponding squark mass-eigenstates, 

9~a = {gi,g 2 }. (22) 

If we neglect intergenerational mixing, that relation is given by 

Q'a = E R S^ 
b 

= ( COS6 l Si i«) (q = t,b), (23) 
y —smVg cos v q J v ' 

where we have singled out the third quark generation, (t, b), as a generical fermion- 
sfermion generation. The rotation matrices in ( fZ3D diagonalize stop and sbottom mass 
matrices: 

M*- ( Ml L +mj + co S 2(3(l-ls 2 w )Ml m t Ml R \ 

Mi ~ \ m t M{ R M? r +m 2 t +l cos 2(3 s 2 w Mf .) ' 

KA 2_ ( Ml L +ml + C os2P(-\ + \s 2 w )Ml m b M b LR \ 

^ m h M b LR M? r +m\-\ cos 2(3 s 2 w M 2 Z , ) ' W 

with 

M£ fl = A t -/xcot/3, M^ = A-/itan/3, (26) 

\i being the SUSY Higgs mass parameter in the superpotential, A tt b are the trilinear soft 
SUSY-breaking parameters and the Mq h R are soft SUSY-breaking masses ||. By SU (2) l- 
gauge invariance we must have Mi L = Mg , whereas Mj R , Mg are in general independent 
parameters. With regard to supersymmetric fermionic partners, from the higgsinos and 
the various gauginos we form the following three sets of two-component Weyl spinors: 

rf = {-iW + ,H+} , Tr = {-tW-,H{} , (27) 

T° a = {-iB^-iWlHim} . (28) 

These states also get mixed up when the neutral Higgs fields acquire nonvanishing VEV's. 
The "ino" mass Lagrangian reads 

C M = - < T + \M\T- >~\< r°|^°|r° > +h.c. , (29) 
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where the charged and neutral gaugino-higgsino mass matrices Ai,Ai° are also well- 
known; in our notation they are given explicitly in Ref.||23||, where we remark the presence 



of the parameter \i introduced above and of the soft SUSY-breaking Majorana masses M 
and M', usually related as M'/M = (5/3) tan 2 (9^. The corresponding four-component 
mass-eigenstate spinors, the so-called charginos and neutralinos, are the following^: 



and 



^ = ( Z a f T J ) = ■ (31) 



]\T* pi 



where the matrices U, V, N are defined through 

U*MV ] = diag{M u M 2 } , N*M Q N ] = diag{M°, ...M 4 } . (32) 

With the help of these matrices, the following interaction vertices appear in the SUSY 
three-body decays of the top quark (after rewriting them in the mass-eigenstate basis for 
all sparticles): 

• fermion-sfermion-(chargino or neutralino): 



C = -gilxT {A + t ai e i P L + AJ ai P R )b 
-gb:xt{A + JPL + AJ a ^P R )t 

Xl U +a fP L + A^e a P R ) t (33) 
xi (A+WPl + AJTe a P R ) b 

+ h.c. 

Here Pl,r — § (1 T 7 5 ) are the chiral projectors, and we have defined the following 
coupling matrices: 



A, 


t 

~ai 


r>{t)*TT* \ r>(t)*TT* 
U la U il ~ A t^2a u i2 


A_ 


t 

~ai 


—Rla AfoVi 2 


A, 


b 

~ai 


K la V il - X b tt 2a V i2 


A. 


b 

~ai 


—R\a ^tUi2 


A, 


(0)* _ 
-act 


Rf; (iV* 2 + Y L t a n6 w N* al ) + R { 2 f V2\ t N* a3 


A. 


(0)t _ 
" act 


R? a *V2\ t N a3 + Rf a *Y t R t^e w N al 


A, 


(0)6 _ 
-act 


42* (K2 - YHanewN*,) - R { 2 fV2X b N* aA 


A. 


(0)6 _ 
- act 


-R ( S*y/2X b N a4 - R®*Y h R %tme w N al , 



(34) 



where Y L and Y/l are respectively the weak hypercharges of the left-handed SU(2) l 
doublet and right-handed singlet fermion-sfermion partners within the chiral su- 
permultiplet. The potentially significant Yukawa couplings are contained in the 



3 We use the notation of Rcf . p3[| , namely, first Latin indices a,b,...=f,2 are reserved for sfermions, 
middle Latin indices i,j,...=l,2 for charginos, and first Greek indices a,/3, ... = 1,2, ...,4 for neutralinos. 
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following ratios with respect to the SU(2)l gauge coupling: 

. _ h t m t _ h b _ m b 

Af = — — — 1= , Ah = — — 1= . (do) 

g V2M w sin(3 g V2M w cos(3 

In the previous formulae we have changed the chargino-neutralino basis from {\l/f , } 
to {xt, Xa} an d m the latter we have introduced the coefficients and e a . The rea- 
son for this change of basis is because in the numerical analysis we use real, instead 
of complex, diagonalization matrices U, V, N for charginos and neutralinos. Thus 
we have to compensate for the minus signs that may appear in the list of mass 
eigenvalues by introducing the e parameters. For instance, for chargino masses we 
set, for each i — 1, 2: 

M i = e i \M i \. (36) 

In this formalism the physical chargino, xt •> is n °t always the the mass-eigenstate 
spinor defined in eq. ([!(]), but 

+ _ f *f if Ci = 1 . . 

X *-\±l^t if Q = -l. (3?) 

This is equivalent to replace 

*r x ~ =(p L + ei p R )^r =Cx t +T ■ (38) 

Indeed, with this proviso the new kinetic Lagrangian for the chargino has always 
the correct sign: 

C K = f + (ip - = xt {ip -\Mi\) xt ■ (39) 

We proceed similarly with neutralinos, for each a — 1, ...,4: 

M a = e a |M a |, (40) 

and 

*° - X° = (Pji + 6 q Pl) • (41) 

Notice, however, that in the real formalism just sketched the physical neutralinos 
are no longer Majorana particles in general (Cf. eq.Q3"T|)), since they satisfy 

X°: = CfJ = e aX ° a ■ (42) 

Once these definitions have been made one has to propagate them carefully over 
all the interaction terms and keep track of the e parametersFI. There are a few 



4 



This procedure is equivalent to the one defined in Refs.|T7], |42fl . We have corrected some sign errors 
and some missing e's in the Feynman rules given in these references. 
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more or less known subtleties related to Majorana particles in connection with the 
e formalism which are worth remembering. Thus e.g. take a generical Lagrangian 
interaction involving a neutral Higgs and two neutralinos: 

C = gH x^5T,sPLX 5 + h.c. 

= ^ x°(e 5 r 75 P L + e 7 q 7 P fl ) X o. t 

Since Xa can be created or destroyed by any of the fermionic field operators ap- 
pearing in the Lagrangian, for a term with fixed indices a and f3 on eq.(fE|) we 
have 

C aP = gH x a (e^ a pP L + e a T^ a P R )x$ 

+gH%(e a r /3a P L + e l3 r* a pP R )x a . 
Therefore, using eq.(f4"2"|) and the usual properties of the charge conjugation matrix, 
we may rearrange the two terms in (|44"D as follows: 

c a p = gH°xi + r> a )P x + e a (T} a + r* af3 )p R ) X J . (45) 

Notice that, for virtual particles, the rule fl38f) just entails the following replacement 
in the numerator of the chargino propagator: 

t+\M i \^>t + e i \M i \=t + M i , (46) 

and similarly for neutralinos, so that one can use real matrices U, V, N together 



with positive or negative mass eigenvalues, since the e's just cancel out[23|. For 



real sparticles in the final state, a similar thing occurs for the square of the various 
amplitudes involved in the top quark decay, but here one has to keep track of the 
e's anyway since they may play a crucial role in the computation of the interference 
terms. In our formulae we shall nevertheless maintain the original complex notation 
so that one can either use complex mixing matrices and forget altogether about e's 
or, alternatively, to keep the e's and consider that all the coupling matrices are real. 

quark-squark-gluino : 

£ = 9 r (A r ) w (4f Pl - R ( £*Pr) Qi + h.c. (47) 



V2" 

where A r are the Gell-Mann matrices. This is just the SUSY-QCD Lagrangian 
written in the squark mass-eigenstate basis. 

squark-squark-Higgs: 

£ = -jfif-ffi* R { : ] c *RSga b H + + h.c. (48) 

where we have introduced the coupling matrix 

/ M^/ sin 2/3 — (m| tan/3 + m| cot/5) — m b (fi + A b tan 0) \ , . 
9ab ~ [ -m t (p + A t cot 0) -m t m b (tan /3 + cot 0) j ' ' ' 
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chargino-neutralino-Higgs: 

C = -gH + xf (sin (3 Q**e a P L + cos (3 Q^Pr) xl + h.c. , (50) 

with 



Qai = VhNca- -^(N a2 + tan6 w N al )V l2 . 



(51) 



Gauge interactions: In our calculation we only need the sparticle interactions with 
the W ± : 



quarks: 
squarks: 



£ = JL frfP L bW+ + h.c. (52) 
v2 



C = iJLBSt B$»W+% # k + h.c. (53) 



charginos and neutralinos: 

£ = 9 X° a Y {C L m e a e % P L + C*P R ) xt W~ + h.c. (54) 



C L m = j= 2 N a3 U* 2 - N a2 U*, 



(55) 



There are of course some additional interactions that can be involved in the SUSY 
three-body decays of the top quark, such as the various quark- quark-higgs vertices, and 
also the vertices involving one gauge boson and two higgses or viceversa. However, they 



are pretty well-known and we shall not quote them here [17 ] . 

4. Numerical analysis 

As mentioned in the preliminary survey of Section 2, among the potentially relevant 
SUSY three-body modes (|i~9"D, only a few result in sizeable corrections to the top quark 
width. The relevant modes are defined to be as those that can give a contribution 

x-pSUSY 

-p (56) 

with respect to the the canonical width @ of the top quark in the SM. For comparison, 
let us recall the size of the one-loop corrections to the canonical decay (0) -these quantum 
corrections being of the same order in perturbation theory as our three-body decays at 
the tree-level. The corresponding SM electroweak corrections lie in the ballpark of 1.5% 
(in the GV-scheme) 0] for the present values of the top-quark mass, whereas the QCD 



corrections are of order —(8 — 9)% [34| and are essentially independent of m t . Our analysis 
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shows that the SUSY three-body decays on eq.(|T9"D that most likely could fulfil condition 
(H) are the modes I, IV, VIII, IX, X. 

As an independent set of SUSY inputs we introduce a similar tuple of parameters as 
in the numerical analyses of Refs. [[□]]- specifically: 

{t^nj3 1 m H +,ii 1 M 1 mg 1 mi v mi v Ml Rl M t LR ) . (57) 

For the full numerical analysis of the ten decays we have produced several hundred 
plots in order to explore all the relevant peculiarities of the 9-dimensional parameter space 
of eq. ( |5T| ) jPJ. In what follows we limit ourselves to report on the main results obtained 



under more simplified assumptions which do not alter the maximal rates expected. In 
particular, throughout all our numerical analysis we use a representative set of inputs 
which is non-critical, i.e. such that our results behave smoothly under variations around 
these values. The simplifying assumptions mentioned above are the following. We shall 
assume that the mixing angle in the stop sector is 9 t = tt/4. In the sbottom sector, 
where mixing effects are more unlikely, we fix M LR = and = = m~ h . With 
these settings the stop mass eigenvalues raj are determined and hence a more restricted 
set of parameters is obtained than the general set (f)7]). In most processes, this suffices. 
Notwithstanding, we shall come back to the general set ( f)% later on for some of the most 
promising decays. 

Apart from the phenomenological limits on the various sparticle masses mentioned 
above, we shall explore the crucial parameter tan/3, eq. (f20|) , within the range 



0.5 <tan/3< 70, (58) 



which is fixed by the perturbative domain of the Yukawa couplings (p5|) . Apart from the 
well established SM inputs |[43| , we adopt the following for our numerical analysis: 

m h = 5 GeV 

m t = 175 GeV , , 

V tb = 0.999 1 } 

a s (m t ) = 0.11 . 



• Decay I: There are two Feynman diagrams contributing to this process (Cf. Fig.l). 
Only two parameters from the general set ((57|) are needed in this case, namely 

(tan P,m H +) . (60) 

Diagram (a) in Fig.l is the conventional diagram, whereas diagram (b) gives the ex- 
tra Higgs contribution, which is determined by the strength of the Yukawa couplings 
given on eq. (|55]) with m T in place of m^. Notice that in this case the charged Higgs 
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need not to be supersymmetric but just the charged member of a general two-Higgs 



doublet model of type II [[17]]. Since m T > m q for any light quark, in the region of 
tan/3 > 1 we do not consider the similar decays in which (z/ T ,r + ) is replaced with 
(u, d) or (c, s). On imposing the condition that the two-body decay t — > H + b is not 
kinematically permissible, we study the three-body decay only for Higgs masses 

m H + > m t — m h . (61) 

Since we aim at a departure from the SM contribution, we define the quantity 

, _ £ft-> bu T T + ) - Tft-> bW + * -> bu T T + ) 

dl ~ T (t^bW+* ^bu T r+) ' [bZ) 

and plot contour lines of 5j in the plane Q60D under the condition flB"I|) . The result is 
displayed in Fig. 9. Notice that V (t — > bv T r + ) on eq.(^) is computed from the sum 
of the two amplitudes in Fig.l, whereas r (t — > bW + * — > bv T r + ) includes the first 
amplitude only. It is seen from Fig. 9 that corrections of a few percent are possible 
at high tan (3. In the low tan (3 < 1 region, the decay under study is inefficient since 
the r-lepton Yukawa coupling becomes very depleted. In spite of the fact that in 
this region of tan/? the alternative decay t —>■ bcs can give contributions of order 
1%, this channel would be much more difficult to separate from the background. 

The dominance of the the semileptonic channel for tan/3 > 0(1) occurs for both 
real or virtual Higgs decays. For real decays, we have 

T(H + -> t + v t ) 1 fm T \ 2 tan 2 /3 



F(H+^cs) 3 \mj {m 2 s /m 2 c ) 2 tan 2 (3 + cot 2 (3 

L / TYl \ ^ / 

-> - — > 10 (for tan (3 > ^m c /m s > 3) . (63) 



3 \m s 

The identification of the r mode could be a matter of measuring a departure from 
the universality prediction between all lepton channels. Fortunately, r-identification 
is possible at Tevatron j|5] ; and the feasibility of tagging the excess of events with 



one isolated r-lepton as compared to events with an additional lepton has been 



substantiated by studies of the LHC collaborations . It has recently been shown 
that it should be fairly easy to discriminate between the H^-daughter r's and the 
if ± -daughter r's by just looking at the opposite states of r polarization resulting 
from the W and H decays; the two polarization states can be distinguished by 
measuring the charged and neutral contributions to the 1-prong r-jet energy (even 
without identifying the individual meson states) 



Decay II: There are four Feynman diagrams contributing to this decay, and they 
are displayed in Fig. 2. In principle, h° in these diagrams could be replaced with any 
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of the neutral MSSM higgses. However, the following relations must hold simulta- 
neously 

m h o < m t -m b - M w , , 

m H + > m t — mfe . 

in order to allow the three-body decay and at the same time to forbid the two-body 
charged Higgs mode (|8|) involved in one of the amplitudes. In the MSSM, the second 
relation implies (at the tree-level) 0] that m^o > 155 GeV, and since the heavy 
CP-even Higgs satisfies m#o > Mz, the relations (U) can only be fulfilled by the 
lightest CP-even Higgs h°. Again the only two parameters involved are those in 
eq. (|60|) , since the mixing angle, a, between the CP-even higgses is not independent 
in the MSSM [[L7[]. Unfortunately, a fully fledged calculation of the diagrams in Fig. 2 
yields a disappointingly small contribution, as seen in Fig. 10. The reason for it is the 
following. The potentially relevant neutral Higgs Yukawa couplings h° 1 1 and h° b b 
in Figs. 2a and 2c are proportional to cos a/ sin/? and sin a/ cos/3, respectively Jl7 



Since m^o is assumed to be large, then a goes to zero or to — n/2 depending on 
whether tan (3 is large or small, respectively. In neither case these Yukawa couplings 
become enhanced with respect to the gauge coupling. The other potentially relevant 
interaction is the charged Higgs vertex H + 1 b in Fig. 2b, which is sensitive to both 



Yukawa couplings on eq. (p5|) and is independent of a [[17]]. Unfortunately, this vertex 
is also rendered ineffective in the present conditions, since its behaviour at large and 
small tan/? is compensated for by the opposite behaviour of the companion vertex 
W + H + h° in the same diagram, which is proportional to cos(a — (3) — > 0. In the 
end there is a balance between the four diagrams giving rise to a maximum near 
tan/? = 1.5, although the resulting contribution with respect to the canonical decay 
is rather meagre, namely 

T(t ^h°bW~ 



making the observation of this process essentially hopeless 



<3xl(T 4 , (65) 



Decay III: The corresponding Feynman diagrams are in Fig.30. Obviously, this 
decay can only proceed within the context of the so-called light gluino scenario 
(9(1) GeV [HI, since for heavy gluinos the usual bound flTBD prevents the 



m 



decay from occurring. With rrig assumed to be negligible, we are left with the 
additional conditions 

(66) 



m~ b < m t — M w 



m~ tl > m t , 



5 We remark that in all our Feynman diagrams indices for virtual sparticles are understood to be 
summed over. 
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where the second relation is to guarantee that the two-body decay (||) does not take 
place. The relevant subset of parameters in the present instance is 

(tan 0,m 5 ,m 6> A4 B ). (67) 

The numerical analysis is presented for a light gluino mass rrig = 1 GeV and fixed 
tan/3 = 1 since small values of tan/3 are preferred. The results are basically the 
same within the light gluino range rrig = (1 — 10) GeV and are rather insensitive 
to allowed values of M t LR below m t , so we have set M l LR = 0. In Fig. 11 we display 
the contribution of the LH sbottom final state. For LH sbottom masses respecting 
the absolute LEP bound fllQ|), we find 

4„^<».., (68) 

1 SM 

which is well below the limit fl5E| ) and therefore is too small. The yield from the RH 
sbottom final state is even smaller (< 10 -5 ) due to the helicity flip at the bottom 
line in Fig. 3a. 



Decay IV:This process is rather complex since it involves four Feynman diagrams 
(see Fig. 4) and the following parameters: 

(tan/3, m H +ilJ>, M, M£ fl ) . (69) 

The simultaneous conditions to be required in order that this decay is possible while 
the two-body modes @-© are forbidden are the following: 

M x o + M x + + m b <m t 
m H + > m t -m b 
M x o + mi > m t 
M + + m~ h > m t . 

Notice that in practice the handling of these relations has to be carried out nu- 
merically since our inputs (|69"D are not given directly in terms of the chargino and 
neutralino masses. Therefore one has to diagonalize the chargino and neutralino 
mass matrices on eq.([29|) and look for regions in the higgsino-gaugino parameter 
space (/i, M) where the relations (|70|) are met. We have checked that the maximum 
contribution of this decay is obtained near the phenomenological boundaries defined 
by the condition (0). In Fig. 12 we plot the ratio 



J IV 



(71) 



as a function of tan (3 for a typical (/i, M) boundary point corresponding to 
60 GeV. We see that at high tan (3 there are significant enhancements of 8iy which 



■xt 
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could boost this quantity up to a value ~ 6 — 8%, i.e. up to the level of the 
conventional QCD corrections to the canonical decay. We have checked that the 
total effect from the next-to-lightest charginos and neutralinos allowed by phase 
space amounts to an additional contribution of 0.6% in the high tan/? region. In 
view of these results, it follows that this three-body mode could be relevant. 



Decay V: Three Feynman diagrams contribute to this decay (see Fig. 5) and the 
parameters involved read 

{Un(3iiM,M,m~ g ,m h Ml R ), (72) 
which are limited by the conditions 



mi < m t — m& — m~ h 

(73) 



nig > m t — mi 



M x o + mi > m t 
M x + + m~ b > m t . 

Since gluinos can either be very light, viz. of 0(1) GeV, or at least as heavy as 
100 GeV, the second condition above enforces a heavy gluino scenario. Within 
our current set of hypotheses we exclude a heavy stop and a light gluino since the 
(degenerate) sbottoms would be heavy, too, and the phase space would be exhausted. 
In Fig. 13 we present a "lattice plot" sampling |fnj of the 6-tuple ([72]), separately 
for LH and RH sbottoms. We see that the ratio 

°v = p , ('4) 

J- SM 

shows the highest cumulative number of points for values 5y < 0.2% whereas higher 
values get only a few number of spots. 



Decay VI: This decay is particularly simple, for it admits a single Feynman diagram 
(Fig. 6a). The input parameters are 

(tan /5,/i,M, m~ b ,m rs ) , (75) 

under the conditions 

2m b <m t -m c 
m- b + M x + > mt, 

where we assume m~ b = m s -. Thus we are led to an scenario of heavy charginos and 
relatively light squarks. In Fig. 14a we display the ratio 

r(t -> b a cS b ) 
ovi = = , (77) 

J- SM 
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for the different combinations of chiral species of squarks. Contributions of order 
Syi — 1% obtain for tan/? < 0.6, i.e. near the lower limit (|58f) . 

Decay VII: This decay is similar to the previous one (see Fig. 6b), but it involves 
r-sleptons in the final state. Thus the inputs are 



bounded as follows 



The relevant ratio 



(tan/3, //, M, m^, ) , (78) 

m~ b + m f < m t 
mr + M+ > m t . ^ ' 



Svn= E&^im , ( 80) 

1 SM 

is analyzed in Fig. 14b for the different combinations of chiral species of staus 
and sbottoms. In the present instance, where sleptons are around, we assume a 
corresponding mass matrix with the same simplified structure as that of the sbottom 
mass matrix. Therefore, the LH and RH staus are degenerate, with a common mass 
m f > 95 GeV, which guarantees rriy > 40 GeV - as required by the data on the 
invisible Z-width || . In these conditions, it is seen from Fig. 14b that Syu — 1% 
is achieved for tan/3 very near the two extreme values of the allowed interval fl5"8|), 
but only at the expense of a rather light chargino M x + ~ 50 GeV, which is still 



permitted if it is almost degenerate with the lightest neutralino |38[ . 



Decay VIII: This decay is related to number IV in that a gluino substitutes for a 
neutralino. The number of Feynman amplitudes stays the same (Cf. Fig. 7) and 
the relevant tuple of independent parameters can be chosen as follows: 

(tan/3,m § ,/i,M,m s ,Mi B ), (81) 

being subdued by the conditions 

M. ' + + rrig < m t — m b 

mi + trig > m t (82) 
m- b + M x + > m t , 

whose meaning should by now be pretty obvious. In contrast to decay V, both 
the light and heavy gluino scenario is permissible in the present instance. In the 
first case (light gluinos), a heavy stop and a heavy chargino could coexist with a 
relatively light sbottom, while in the second case (heavy gluinos) a light stop is 
tolerated at the expense of heavy sbottoms. The numerical analysis of the quantity 

6v,u= T(t Z b ~ gXt) , (83) 

J- SM 
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shows that in the heavy gluino scenario the contributions are well below 1% and we 
do not show explicit plots. In contrast, for light gluinos we see in Fig. 15 that Syin 
could border on values of order 1% for any value of tan /3, and it could even reach 
4% for sufficiently small or big values of tan/? within the interval ([58]). This result 
is quite remarkable, since it is about 60% of the QCD corrections to the canonical 
decay. 



Decay IX: This is a very interesting three-body decay to deal with, though its 
analysis is technically quite demanding for it involves four Feynman amplitudes 
(see Fig. 8). In the Appendix at the end of the paper we present the full squared 
amplitude corresponding to this particular decay. The parameter space includes as 
many as nine parameters, namely 

(tan p, m H + , //, M, m s , m~ bi 2 , m~ tx , M{ R ) . (84) 

In this case, and due to the potentiality of this mode, we relax the assumption on the 
stop mixing angle being 9 t = 7r/4 and introduce the lightest stop mass as a new 
input. Furthermore, we also abandon the restriction M b LR = and the assumption of 
degenerate sbottom masses by making allowance for a free input value of the sbottom 
mass eigenvalues mg with the proviso that 6b = 7r/4. In contrast to the previous 
decays, we shall not impose conditions blocking the possibility that the intermediate 
two-body states in Fig. 8 can be real two-body decays, except for the gluino mode 
t tg which, if allowed, would be overwhelming in most of the parameter space. 
As for the remaining two-body modes, it turns out that the present three-body 
decay can be competitive with them in certain regions of parameter space. 

The numerical analysis of decay IX confirms the expected fact that it can be relevant 
only if the two-body channel H + — > tb is kinematically forbidden, otherwise the 
Higgs decay width becomes too large and it has a dramatic suppression effect on the 
partial width of our three-body decay. This can be seen in Fig. 16a, corresponding 
to mi = 60 GeV, m~ b = 100 GeV, rrig = 130 GeV and a large value of tan/3 of 
order m t /m b ~ 36. We have also fixed fi = M = 100, 250 GeV and M b LR = 0. We 
see indeed that, for m#+ approaching mi + m- b from below, the decay IX can give a 
contribution to the ratio 

f n = ^, (85) 

J- SM 

which ranges from a few percent to 100% and above. The contribution from 6 2 is 
basically the same. Clearly, this behaviour is unmatched so far by any of the previous 
decays. Therefore, these non-standard effects are to be seriously considered in any 
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consistent treatment of the total top quark width. We postpone for Section 5 the 
discussion of the possible signatures. 

Since in the conditions under study, the two-body mode t — > H + b is still available, 
we also plot in Fig. 16a the ratio 

E<^>. (86) 

J- SM 

In this way we see that, depending on the value of /x, there may be a sizeable range 
of Higgs masses where the decay IX is not negligible as compared to the decay 
t — > H + b. Although these results have been obtained by assuming that the two 
sbottoms are degenerate, we may also use the general set of inputs ([84] ) to enter 



different values for the sbottom masses. For example, take = 100 GeV and 
mi 2 = 150 GeV, and set \x = 100 GeV without altering the rest; we then obtain 
a rate 5jx = 33%. In all these cases the soft SUSY-breaking parameter A b is 
necessarily large to avoid conflict with the phenomenological mass bounds. In fact, 
if in the previous example we would trade the input parameter \x for Ab, we would 
find that we cannot demand a small input value Aj, ~ since it would entail a too 
small a value of /j below 50 GeV, which is ruled out by the chargino mass bound 



Since the general set (|84]) is found not to change the order of magnitude of the rate 
of our decay, it will simplify the analysis to restore the original set of inputs 

(tan/3, m H +, /i, M,nig,m~ b , Ml R ) . (87) 

where 6 t = 7r/4, M b LR = together with the assumption of degenerate sbottom 
masses. In Figs. 16b and 16c we study the evolution of 8jx versus tan/3 for the 
two regimes where the Higgs decay H + — > tb is closed and open, respectively. In 
the closed case the slope is very high in the large tan/3 region. As a consequence, 
the yield from the Higgs mediated diagram becomes so overwhelming at large tan j3 
that the rate of the decay IX overtakes that of the canonical decay as soon as 
tan (3 > 50,!. Notice that the long flat region in the intermediate interval 1 < tan (3 < 
20 is sustained by the gluino mediated diagram in Fig. 8, though in this plateau 
the contribution stays small due to the large mass assumed for the gluino, which 
prevents the t —*■ tg decay from occurring. Whereas the dependence on the soft 
SUSY-breaking parameter M is very mild, in Fig. 16d we exhibit the dramatic 
dependence on the parameter /j for typical values of the other parameters. As 
already mentioned, the region \i < 60 GeV (which is unfavoured by our decay) is 
excluded by the present bounds on chargino masses. Notice also that when the real 
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Higgs decay into tb is viable the ratio Six lessens considerably, viz. down a few 
percent level and mostly at a few per mil or below. 

Intriguingly enough, there are regions of parameter space where the decay IX can 
be competitive with the main two-body modes available. In Figs. 17a and 17b we 
plot the two-body decay rates 

r(t -> h x ° a ) ^ r(t -> kxt) j (88) 

as a function of tan (3 for two values of /i and for fixed M + = 60 GeV. For the 
squark masses we have taken mg = 100 GeV, = 60 GeV. In the first decay 
we have included separated plots for the first three neutralinos xi, X2, X% since they 
are also allowed by phase space. Although the last two are heavier than xi they are 
worth including, the reason being that they are more higgsino-like and therefore have 



larger Yukawa couplings of the form (35) resulting in a competitive contribution. 
For the second decay in (88), the chargino xt nas n °t been included since it is too 
heavy. It is seen that the rates can be high for small tan j3, and in this region they 
are more efficient than Six- However, for tan/5 > 30, the three-body decay width 
Six is n °t om y competitive but it can even surpass the rate of the previous two-body 
decays, especially for large enough \i where the partial width of the former increases 
substantially whereas the two-body partial widths (B8l) decrease. 



• Decay X: This decay is similar to the previous one. The Feynman diagrams are 
similar to those in Fig. 8 by just replacing t — > v T and b — > f and forgetting about 
gluino and neutralino mediated amplitudes. The parameter space for this decay 
reads 

(tan Tfiu+ , ^t, iW^^j, Tfiy T ) . (89) 
An obvious restriction to be fulfilled is 

rriy < m t — nib — m? ■ (90) 

For simplicity we shall assume that the two f-sleptons are degenerate in mass and 
that M T LR = 0. As in the decay IX, we expect that the ratio 

J- SM 

will be most sizeable if the two-body Higgs decays H + — > z> T f + and H + — > tb are 
kinematically forbidden. If this is so, for m#+ approaching + from below, 
the decay X can furnish a contribution to the ratio Sx which ranges between the few 
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percent to the several ten percent, depending on the value of the other parameters, 
specially of tan/3 and /i. This is seen in Fig. 18a, where the relevant inputs are 

(tan/^mpj = (40,100 GeV,50 GeV). (92) 

Cases fi and f 2 are almost indistinguishable since tan/3 is high in Fig. 18a. In 
Figs. 18b and 18c we plot the evolution of 5x on tan/3 and on fi, respectively, for 
fixed values of the other parameters. Notice that on increasing tan/3 the Higgs 
coupling to stau and the corresponding sneutrino increases but at the same time 
the stau mass also increases. For = 50 GeV, this mass saturates at a value 
of ruf ~ 92 GeV for tan/3 ~ 10. Therefore, for larger values of tan/3, the ratio 
5x steadily grows and it can surpass 50% for tan/3 > 55. These results can still be 
improved substantially for higher values of (Cf. Fig. 18c). 



5. Discussion and conclusions 

To summarize, we see that there are some three body decays of the top quark within the 
MSSM whose contribution can be a fraction of the standard width Tsm — r(t — > W + b) of 
order or above 1%. The latter reference value is approximately the size of the conventional 
electroweak quantum corrections to Tsm For the sake of comparison, let us point out 



that the Higgs effects on t — > W + b within the MSSM are about one order of magnitude 
smaller than our reference value, i.e. they are of order 0.1%|[22|. Many of our three-body 
decays could give a contribution above this small number. In general, however, we find 
that only a few number of these three-body decays do contribute, within the same order 
of perturbation theory, as much as the full (electroweak plus QCD) quantum effects on 
t — > W + b in the SM, which are of order —7% ||34|| . Nonetheless, if we consistently add up 
the yield from some of these decays, we are able to find regions of the parameter space 
where the resultant pay-off could counterbalance the one-loop quantum effects. This 
could be the case e.g. if we add the contributions from decays I and IV within the heavy 
gluino scenario, or decays I and VIII within the light gluino scenario. In either of these 
situations the one-loop corrections would appear as "missing effects" in a measurement 
of the total width of the top quark, T t . This feature is remarkable, since it is compatible 
with an scenario in which all of the SUSY two-body decays of the top quark are blocked 
up. In such a case the leading SUSY signature in inclusive top physics could come from 
the combined effect of the available three-body decays. 

Finally, there are two cases, viz. decay IX and decay X, which could be very important 
in certain regions of parameter space, specifically in a region where both tan/3 and \fi\ are 
large enough, tan/3 > mt/mb — 36 and \y\ > 100 GeV, provided that the Higgs mass is 
below the z/ T f + and t b thresholds. In this case, the last two decays could contribute to the 
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total top quark width a fraction which, if one does not stretch too much the values of the 
relevant parameters, it can still range between the few percent to the several ten percent. 
The large size of these contributions are the tree- level counterpart, within the same order 
of perturbation theory, of the large one-loop quantum effects on the two-body decay 
(H) P^, Therefore, decays IX and X could possibly be searched as exclusive decays 
since they may give rise to well defined and rather exotic signatures. These signatures 
have to be compared with the canonical signature 

t - bW + - b(ff') , (93) 

which consists either of 3 jets (one of them a b-jet) or a b-jet + a positively charged 
lepton ( l + )+ missing energy- momentum ( In Table 1 we display the alternative 
decay signatures associated to decay IX. In all cases we have obviated an additional b-jet 
which is in common with the canonical signature. Therefore, what we show is the specific 
signature associated to the virtual decay H + — > t&[|33|. In some of the signatures in 
Table 1 (cases (a)-(c) ) we have admitted of the exotic possibility to having FCNC decays 



mediated by neutralinos [21|. These could enter the game if the chargino decay mode of 
the stop, 

h - bxt , (94) 

is kinematically forbidden. If this mode is allowed, it would be dominant and the main 
signatures would be (d)-(f) in Table 1. As for the sbottom decay, the charged mode 
b — > cxt is CKM-suppressed whereas b — > bxi is always available (if a light neutralino 
exists). Let us point out that if the former case applies it would lead to two final leptons 
with the same electric charge (see case (f)), which is definitely a non-canonical signature. 
In general, we should expect that the leading signatures in Table 1 are (a)-(b) or (d)-(e), 
depending on whether the channel (0) is closed or open, respectively. 

As far as decay X is concerned, a similar discussion ensues. The main signatures are 
detailed in Table 2. It could well happen that the charged sneutrino decay 

Vr^T'xt, (95) 



is phase space obstructed, similarly to the previous case fl94|) for the stop. However, an 
important difference here is that the neutral decay 

Vr^VrXl, (96) 

is perfectly possible and most likely it is the leading decay mode of the sneutrino, if Xi is 

the lightest SUSY particle. Another difference with respect to decay IX is that whereas the 

6 Notice that some of the Yukawa couplings ( |35| ) and triple scalar couplings (^) can be comparable 
to the strong gauge coupling (|47]). Therefore, the tree- level corrections to decay IX can be of the same 
order as the SUSY-QCD p6f and electroweak SUSY corrections [p8| to the two-body decay (ph. 
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charged decay of the sbottom is CKM-suppressed, the charged f-decay f + — > u T Xi is n °t- 
We conclude that the leading signatures are (a)-(c) in Table 2, where the characteristic 
trait is the proliferation of jets, lepton prongs and missing energy. However, if the charged 
decay mode (|95| ) would be available, it could be rather interesting, for the corresponding 
signatures (d)-(f) in Table 2 show the distinctive presence of a r-lepton with the "wrong" 
sign, i.e a r~, in contradistinction to the canonical decay fl93|) where one expects to 
find a r + . Notice that in the conditions under study the signatures on Tables 1 and 2 
cannot come from the decay of a real H + ; in fact, although real charged Higgs bosons 
can be produced in these conditions, their leading decay would be H + — > t + v t . With 
enough statistics, the latter signature would lead to a violation of lepton universality and 
it has a chance to be isolated in the future (Cf. Section 4); however, with an insufficient 
statistics it could easily be confused with a standard signature (|93|). Quite in contrast, the 
exotic signatures displayed in Tables 1 and 2 are highly non-standard and may play an 
important role in the top quark phenomenology beyond the SM, especially in the analysis 
of the single top-quark production processes [[!(]] , which are foreseen to play a decisive role 
in the near future at Tevatron and at the LHC 
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Appendix A 

In this appendix we write out the full analytical expression for the squared amplitude 
of one single decay, since the complete formulae for all the processes that have been 
analyzed is extremely lengthy. We have chosen the decay IX, which is one of the most 
relevant candidates in the list of SUSY three-body modes presented on eq. (|19D . Using 
the Lagrangian interactions and coupling matrices defined on eqs. (147| ) - (149|) , 
(|52]) -(|54|), along with the well-known SUSY Higgs-fermion interactions ||17|| , one computes 
the Lorentz invariant amplitudes Tj (i = 1, ...,4) of the decay t — > btbb a (a,b = 1,2) in 
terms of the four-momenta specified in Fig. 8a. The total squared amplitude summed 
over all final spin, colour and sflavour indices and averaged over the initial states can be 
casted in the following way: 

E<m 2 > = HQu + E {Qa + Qh) (a- 1 ) 

i i<j 

where 

= E(^t) . (A.2) 
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The explicit structure of the various Qij is rather cumbersome. We first do the following 
definitions^ 

6\y = Mw^w 
= m H +T H + 

^ab — 9cd tida K cb — 9LL ti Xa K lb + gRR K 2a K 2 b + 9LR ttla K lb + 9RL ti\a K 2b 

Sum = J2 a ,ij,s,t X ijKt = 'tt 
and derive the following results: 



Qo 



3g*R^ 2 R[f 

•i la lb 



' - 2 

(kqi (pip 3 - p 2 p 3 ) 



2 (pxk - p 2 k) (pxp 3 - p 2 p 3 ) 



- [m~ 2 + mi 2 - 2p 1 p 2 ) p 3 k 



mw 
f rrir 2 — mr' 



mw 



(mi 2 - m~ b 2 ) p 3 k + (pxk - p 2 k) p 3 qx) 



(2 kqi p 3 qi ~ Pzk 



Q 



bb 



3.9 4 G ab 2 

im w 4 [ (qi 2 -m H+ 2 



( 2 



HT J ■ t H+ J 

m 2 m 2 + p 3 k [m 2 cot [3 2 + m 2 tan /5 2 )) 



= ,< T\v hm b m § 2 m t R^R^Rf b ) R^ b ) 
+2m b m t q 2 2 R^RSRf b ) R% 
-Zm~ g m t p 3 q 2 \ K la K lb K 2b + K 2a K lb K 2b I 

fcg 2 m 5 I _K la .ft 2a /% + K la K 2a K 2b J 
+ (2 kq 2 p 3 q 2 - p 3 k q 2 2 ) (it!?/ it*<f + R 2 f R® 2 ) ) 

7 Here, for convenience, we express our results in real, instead of complex, matrix notation following 
the rules stated in Section 3. 
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a, 2 m,,m (a {0)b A (0) * A ^ b A (0) * 4- A {0)b A {0)t A ^ b A {0)t "l 
q 2 m b m t yA +a/3 A +b/3 A ~aa A ~ba + A - a p A -bj3 A +aa A +ba ) 

+ (2 * ft ^ - m g 2 2 ) K$ 6 < )4 a + t + A_i / a4 / 

+^2 m t (A$» A + {f A J b f + A_g» A + £><) 

4-^ 777,777 n (0)6 4 (0)t /l (°) 6 4 (0) '+/| { ° )b A {0)t A ^ b A {0)t \ 

-\-Kq 2 mbm x a [A +a/3 A +b/3 n -aa n +ba + A - a f3 n -bf3 n +aa n -ba J 

+Pm m tm X {) f3 [ A -af3 A +bf3 A ~aa A ~ba + A +af3 A ~b/3 A +aa A +ba ) 

4-^777,7770 (0)fe /4 {0)t A ^ b A {0)t + A {0)b A {0)t A A {0)t \ 
+Kq 2 m b m x o /3 \ A- aj3 A +b/} A +aa A +ba + A +a/3 A_ b/3 A_ aa A_ ba ) 

4777, 777 777 n 777 n (A ^ A {0)t A ^ b A (0) * 4 A A {0)< A ^ b A {0) ^ 

-rrn b m t m x o a m x o [) [A_ a/3 A +b/3 A +aa A_ ba -+- A +af3 ^_ b/3 /i_ aQi /i +6a j 

4-77,^- 777 n 777 n (a {0)b A (0)t A ^ b A {0)t 4 A (0)b A {0)t A A (0) *^ 
+p 3 Km x a m x o /3 \A- a p A +bf) A- aa A +ba + A +a(5 A^ b/3 A +aa A_ ba jj 



Qab + Qab — 



2e H+ e w )R^R$G ab 



3g 4 (2(gi 2 -m g+ 2 ) (qi 2 -m n 

4m w 2 (^( qi 2 -m H+ 2 ) 2 +e h+ 2 ^J ((q 1 2 -m w 2 ) 2 +e w 2 ) 

mi 2 cot (5 [P2P3 - PlP3 + 



+m b 2 tan (3 ( p 2 k — p x k + 




n n j _ g 2 9s 2 (qi 2 -m w 2 ) SumR ( V R$ 

6l(gi 2 -mi4/ 2 ) +e w J (g 2 2 -mg 2 J 



&?2 ( P 2 P3 - PlPZ + "^J 2 ' P3?l ) + (p2k~ Plk + "\„" S k( ll ) 



-p 3 k \ p2q2 -Piq2 + mi mv ^~ h j 
-m b m- g (p 2 k -pik+ mi ~™> h kq^j B$ 

( m- 2 -m- 2 \ D (6) D (t) 

-m- g m t lp 2 P3 ~ PlP3 + mw 2 b PsQl J R\a R2b 

+m b m t \p2q2 ~Piq2 + 9i«2 ) 



^ad + V ac( - 2 ((q 1 2 -mw 2 ) 2 +6w 2 )(92 2 -m x Q 2 ) 

-Psk (p 2 q2 - P^2 + ^£ qi q 2 )) A + [f 
+m b mi {P2q2 - P^ + ftft) A_[f A_^ 

+m b m x% (p2k- Pl k + ^£k qi ) A_ a fA + ^ 
+m t m x o a l-pm+P2P3 + ^p^P 3 qi) A + ^ b 



27 



,f _ g 2 9s 2 (qi 2 -m H+ 2 ) SumG ab 

6m w 2 ^(<ji 2 -m H+ 2 ) 2 +e H+ 2 1 (<?2 2 -mg 2 ) 

m-g p 3 k (—m b tan (3 42 42 — Tn t cot (3 42 4t ) 
+m b m g m t (-m t cot (3 42 4? - tan /3 42 4?) 
+m{, kq 2 (m b tan /? 42 4& + m * c °t /3 42 4b ) 
+m 4 p 3 g 2 (rr^ cot (3 42 #i& + m & tan /3 42 41 ) ) 



( 



m b kq 2 (m b tan (3 A + ^ b A + ( b f + m t cot 
+m t p 3 q 2 (mt cot /3 A+g 6 A+ ( b f + m b tan (3 A. ^ A„ (?) 
+m x o Q M (m 6 tan /? AJ% b A+E* + m t cot (3 A + ^ b AJ b f) 
+m b m t m x o Q (m t cot /? A +b °^ + m b tan /3 A + ^ b AJ b f) ) 



g 2 g s 2 Sum / 

12(g2 2 -mg 2 ) (q2 2 -m x o ^ 2 ) 

-m g m t p 3 q 2 [A- af3 A_ bf3 R 2a R lb + A +af3 A +bp R la R 2b j 

m m m m , f 4 ( ° )6 4 (0) * R^ R W 4- 4 4 (0) * R (6) R ( *^ 

-m h m g m t m x o [A +aj3 A_ bp R 2a R lb + A_ af3 A +bf3 R la R 2b j 

mm n h ( 4 (0)fe 4 (0)t R (6) R W 4- 4 (0)fe 4 (0) * R (6) R ( *^ 
-m g m x o p 3 k [A_ af3 A +b/3 R 2a R lb + A +a(j A_ b(3 R la R 2b j 

m m ha ( A ^ b 4 (°)* R^ E?(*) I 4 4 (°)* r( 6 ) d(«)\ 

-m b m g kq 2 [A +a/3 A +bf3 K 2a K lb +A„ a(] A_ b/3 K la K 2b J 

Mm m a 2 ( A 4 (°)* R^ !?(*) _j_ 4 4 (°)* R^ 
+m b m t q 2 [A- a p A_ bp K la K lb + A +a/3 A +bf3 K 2a K 2b J 

Mm m „ n n ( A ^ b A (°)* f?( 6 ) R^ _L 4 4 (°)* R^ R^ 

+m t ,m x o f3 p 3 q 2 [A +af3 A_ bf3 R la R lb + A- af3 A +bf3 R 2a R 2b ) 

Mm m , ha ( A ^ A (°)* /?(*) /?W _L 4 4 Z?C*)S 

+m b m x o f3 kq 2 {A_ a(j A +b(3 R la R lb + A +a/3 A_ bf} R 2a R 2b ) 

M(9hn n n n b n 2 \ ( A ^ b A ^ R^ pW _L 4 4 (°)* r( 6 ) rW^ 

+ {2kq 2 p 3 q 2 - p 3 kq 2 ) [A +a/3 A +bf] R la R lb + A_ a/3 A_ bf3 R 2a R 2b j j 



The decay rate is obtained from the formulae above by converting all scalar products 
PiPj = pi.pj involved in these equations in terms of Mandelstam invariants Si, S 2 , S 3 and 
then by performing integration over them in the standard manner j|8] . Care has to be 
exercised in the numerical integrations near the poles. 
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Figure Captions 

Fig.l Feynman diagrams for the decay I. Diagram (a) is the SM amplitude, and 
diagram (b) is the extra charged Higgs contribution in the MSSM. 

Fig. 2 Feynman diagrams for the decay II. 

Fig. 3 Feynman diagrams for the decay III. 

Fig. 4 Feynman diagrams for the decay IV. 

Fig. 5 Feynman diagrams for the decay V. 

Fig. 6 Feynman diagrams for the decays VI (a) and VII (b). 

Fig. 7 Feynman diagrams for the decay VIII. 

Fig. 8 Feynman diagrams for the decay IX. Those for the decay X are obtained 
from (a) and (b) after replacing t — > u T and b — > f. 

Fig. 9 Isolines of 5j in the tan (3-m H + plane. 

Fig. 10 5u as a function of tan/5 for two values of m H +. 

Fig. 11 5 in as a function of m~ b for a b L final state, with tan (3 = 1 and m~ g — \ GeV. 

Fig. 12 Maximum value of 5iv as a function of tan/?, for m~ h = 125 GeV, mi = 
155 GeV, m H+ = 175 GeV and m v+ = 60 GeV. 

Fig. 13 Lattice plot sampling of S v for b L and b R final states as a function of 
m~ b , with tan/3 = 40, m g = 140 GeV, -5M Z < /i < 5M Z , < M < 5M Z , and 
125 GeV < M\ R < 230 GeV. In this case we have used m x + > 50 GeV. 

Fig. 14 (a) Maximum value of Syn as a function of tan/3 for m x + = 135 GeV and 
m~ h = m~ a = 50 GeV; (b) As before, but for Sym and m f = 95 GeV. 

Fig. 15 Maximum value of Svin for fixed m g = 1 GeV, m t ~ = 190 GeV and m~ h = 
150 GeV. The maximum is attained at points on the (/z, M)-plane corresponding 
to the phenomenological limit m + = 60 GeV. 

Fig. 16 (a) Six as a function of m H + for tan/3 = 36, m t ~ = 60 GeV, m~ b = 
100 GeV, m g = 130 GeV and /i = 100 GeV, 250 GeV. Also plotted is the ratio 
T(t — > bH + ) /r sm] (b) Six versus tan/3 for /i = 100 GeV and a charged Higgs mass 
value below the tb threshold; (c) As in (b) but for a charged Higgs mass above 
the tb threshold; (d) Six as a function of //, for the same inputs as in (b) and 
tan/3 = 36. 
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• Fig. 17 (a) T(t — > tXa)/^SM for the first three lightest neutralinos as a function of 
tan/5 for m % + = 60 GeV, m~ t = 60 GeV and fi = 100 GeV, 250 GeV; (b) As before, 
but for the decay t — > & a Xi~ and = 100 GeV. 

• Fig. 18 (a) 5 X as a function of for tan/3 = 36, = 50 GeV and \x = 

100 GeV; (b) 5 X versus tan/3 for m H + = 140 GeV; (c) S x as a function of fx for 
the same inputs as before. 



Typical Decays 



Signal 



(a) 



ex? 



h,2 -> 6x? 



2 j's + jj 

(more ^, less & activity) 



(b) 



6i, 



b{z^)xi^b{ff) x \ 



or 



6(ft°xS)-5(66 or r+r-) X ? 



4 j's + # 2 j's +/$; 2 j's + /+ r + j}. 

(Z° or /i° emission; 

more 6 activity if ft. or Z° — > bb; 

less 6 activity if h° or Z° — » // or vv) 



(c) i I 1 ~* CX? 

1 &2 - cx+ - c(//'x?) 



4 j's 2 j's + 1+ +j>. 
(isolated lepton l + ) 



, . , ii - &xi" - H/f x?) 

' 1 &i,2-6x? 



4 j's 2 j's + 1+ +j>. 
(more more 6 activity) 



(e) 



ii - &Xi~ - K/fx?) 

6(ZW)xS) - 6(//)xS 
&1.2 -^bx^ _or 

b{h\\)^b{bb or t+t-)x? 



6 j's + jj; 4 j's + l+l~ + i>; 4 j's + /+ + jj; 

4 j's + jj; 2 j's + r+ + # 2 j's + /+ + jj. 

(Z° or /i° emission; 

more 6 activity if ft or Z* — » 66; 

Trilcpton) 



(f) 



ii - &xi" - K/fx?) 

61,2 - cx+ - c(/fx?) 



6 j's + # 4 j's + Z+ + i>\ 2 j's + 1+ l'+ + f>. 
(dilepton of same sign) 



Table 1: Detection signatures for the decay IX. The decay chains (a)-(c) are possible only 



ifm fl < m^ i - i 



and m x o < m^. The decay chains (d)-(f) are possible only if > 



xi 



Here ^,j,/" ± ,Z ( *) and / denote missing energy-momentum, jet, isolated charged lepton, Z° 
real or virtual, and (q, ^,1^), respectively. 
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Typical Decays 



Signal 



(a) \ Ut ~" VrX " t+ +j> 

1 Tr2 -> T+xl (more /}) 



(b) 



1,2 T X5 



r+(Z(*)) X ? 
or 

r+(/i°X?) - 



r + (//)Xi 
"(66 or T+- 



)X? 



2 j's + r+ + # t+ + Z+r +f 

(Z° or h emission; 

more 6 activity if h or Z° — > 66; 

less 6 activity if ft or Z° — > // or i/z/) 




2 j's + /}; /+ + j>. 
(isolated lepton; 
less 6 activity) 

2 j's + r+ t~ + j)\ l + + r+ r~ + j*. 
(more /}) 



(e) 



T Xl" 



■(//'x?^ 



-+ 



(^ W )x?) 



r+(//)x? 



or 

r+(^°X?) 



t+(66 or r+r-)x? 



4 j's + r+ t" + i>\ 2 j's + r+ r" + Z+ Z~ + # 

2 j's + t+t~ + /+ + # 

r+ t- + /+ + l'+ I'' + j). 

(Z° or h° emission; 

more 6 activity if h° or Z* — > 66; 

trilepton) 



(f) 



"(/fx?) 

*rtf - M//'X?) 



T Xl" 



4 j's + t- + # 2 j's + t~ + 1+ + i>: 

t- + 1+ 1'+ + ■/>. 

("wrong" -sign single r lepton; 
same sign dilepton) 



Table 2: Detection signatures for the decay X. The decay chains (a)-(c) are possible only 
if my > m x o. The decay chains (d)-(f) are possible only if rriy > rn x +. Here ^ ) },l ± ,Z^ 
and / denote missing energy-momentum, jet, isolated charged lepton, Z° real or virtual, 
and (q, l ± ,i'), respectively. 
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